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(54) Air-fuel ratio control system for multi-cylinder internal combustion engines 



(57) An air-fuel ratio control system for a multi-cylin- 
der engine is provided. An air-fuel ratio sensor is 
arranged in an exhaust system of the engine for detect- 
ing an air-fuel ratio of a mixture supplied to the engine 
and for generating an output indicative of the air-fuel 
ratio of the mixture. An adaptive controller determines 
an amount of fuel to be supplied to the engine with a first 
predetermined repetition period in a manner such that 
the output from the air-fuel ratio sensor becomes equal 
to a desired value. An adaptive parameter-adjusting 
mechanism adjusts adaptive parameters used by the 
adaptive controller. In the adaptive parameter-adjusting 



mechanism, the adaptive parameters are calculated 
with a second predetermined repetition period longer 
than the first predetermined repetition period, and out- 
put data indicative of results of the calculation is gener- 
ated. Further, the output data indicative of results of the 
calculation is smoothed and output data indicative of the 
smoothed data is generated with a repetition period at 
least equal to the first predetermined repetition period. 
The adaptive controller uses the smoothed data as val- 
ues of the adaptive parameters. 
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BACKGROUND QF THE INVENTION 
Field of the Invention 

[0001] This invention relates to an air-fuel ratio control system for multi-cylinder internal combustion engines, which 
controls an amount of fuel to be supplied to the engine by means of feectoack control based on an adaptive control the- 
ory. 

Prior Art 



[0002] Conventionally, an air-fuel ratio control system for internal combustion engines is known, e.g. from Japanese 
Laid-Open Patent Publication (Kokai) No. 8-232733, in which an adaptive parameter-adjusting mechanism calculates 

15 adaptive parameters and an adaptive controller carries out feedback control of the air-fuel ratio of a mixture supplied to 
the engine to a desired air-fuel ratio by using the adaptive parameters. In the known air-fuel ratio control system, an air- 
fuel ratio sensor arranged in the exhaust system of the engine detects the air-fuel ratio of the mixture and supplies a 
signal indicative of the detected air-fuel ratio to the adaptive controller, which in turn determines an amount of fuel to be 
supplied to the engine such that the detected air-fuel ratio becomes equal to a desired air-fuel ratio. 

20 [0003] In the air-fuel ratio control system disclosed in the above publication, values of the adaptive parameters calcu- 
lated in the past by the adaptive parameter-adjusting mechanism are averaged, and the adaptive parameter values thus 
averaged are used by the adaptive controller, so as to prevent the adaptive control from being adversely influenced by 
a particular cylinder when the adaptive control is carried out in a manner synchronous to the combustion cycle of the 
engine. 

25 [0004] Further, if the calculation of a plurality of adaptive parameters is carried out with the same repetition period as 
a repetition period of calculation of the amount of fuel to be supplied to the engine, e.g. a repetition period of generation 
of TDC pulses, the amount of calculation increases to increase toad on a CPU which calculates the adaptive parame- 
ters. To avoid this inconvenience, the repetition period T8 of calculation of the adaptive parameters is set longer than 
the repetition period TF of calculation of the amount of fuel to be supplied to the engine, e.g. to a period equal to TF x 

30 the number of cylinders, to thereby reduce the amount of calculation of the adaptive parameters to reduce the load on 
the CPU. 

[0005] In the conventional air-fuel ratio control system, however, when the averaged adaptive parameter values are 
calculated by the adaptive parameter-adjusting mechanism with the repetition period T8 of calculation set longer than 
the repetition period TF with which the adaptive controller calculates the amount of fuel supplied to the engine, a plu- 

35 rality of values of each adaptive parameter calculated at the repetition period TO at different times are averaged and the 
averaged parameter values are used by the adaptive controller. This brings about the following problem: 
[0006] The adaptive parameter values used by the adaptive controller each contain a component of a particular fre- 
quency corresponding to theteciprocal of the repetition period Te and harmonic components of the same, whereas the 
output from the air-fuel ratio sensor does not contain components corresponding, especially, to the harmonic compo- 

40 nerrts due to the low-pass characteristic of the sensor. This increases an identif ication error in the adaptive parameters, 
resulting in divergence of the adaptive control. 

[0007] More specifically, the adaptive parameter-adjusting mechanism is capable of carrying out accurate identifica- 
tion of the adaptive parameters on condition that frequency components of output data (= the adaptive parameters) 
from the adjusting mechanism are fed back as input data to the adjusting mechanism without being lost after being 
45 transformed through a controlled variable (= the amount of fuel to be supplied to the engine), a plant to be controlled (= 
the multi-cylinder engine), and a plant output sensor (= the air-fuel ratio sensor). However, in actuality the components 
corresponding to the harmonic components are lost due to the low-pass characteristic of the air-fuel ratio sensor and 
not fed back to the adaptive parameter-adjusting mechanism. This causes the above-mentioned increase in the identi- 
fication error. 

so [0008] Further, since the repetition period T6 of calculation of the adaptive parameters is set longer than the repetition 
period TF of calculation of the amount of fuel to be supplied to the engine, the calculated adaptive parameter values 
cannot be quickly changed in response to changes in operating conditions of the engine without delay, which hinders 
the amount of fuel to be supplied to the engine from being quickly converged to a desired value. This inconvenience 
becomes more marked through averaging of adaptive parameters since the averaged values of the adaptive parame- 

55 ters are applied to the calculation of the amount of fuel with the repetition period Te. 
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SUMMARY OF TH £1NVEN^P 



[0009] It is an object of the present Invention to provide an air-fuel ratio control system for a multi-cylinder internal 
combustion engine, which is capable of continuing stable adaptive control even when the adaptive parameters are cal- 
culated with a repetition period longer than the repetition period of calculation of an amount of fuel supplied to the 
engine. 

[0010] To attain the above object, the invention provides an air-fuel ratio control system for a multi-cylinder engine 
having a plurality of cylinders and an exhaust system connected to the cylinders, including air-fuel ratio-detecting 
means arranged in the exhaust system for detecting an air-fuel ratio of a mixture supplied to the engine and fey gener- 
ating an output indicative of the air-fuel ratio of the mixture, an adaptive controller for determining an amount of fuel to 
be supplied to the engine with a first predetermined repetition period in a manner such that the output from the air-fuel 
ratio-detecting means becomes equal to a desired value, and adaptive parameter-adjusting means for adjusting adap- 
tive parameters used by the adaptive controller. 

[001 1 1 The air-fuel control system according to the invention is characterized in that the adaptive parameter-adjusting 
means comprises: 

adaptive parameter-calculating means for calculating the adaptive parameters with a second predetermined repe- 
tition period longer than the first predetermined repetition period and for generating output data indicative of results 
of the calculation; and 

smoothing means for smoothing the output data from the adaptive parameter-calculating means and for generating 
output data indicative of the smoothed data with a repetition period at least equal to the first predetermined repeti- 
tion period; 

wherein the adaptive controller uses the smoothed data generated from the smoothing means as values of the 
adaptive parameters. 

[001 2] Preferably, the smoothing means stores the output data from the adaptive parameter-calculating means with 
the first predetermined repetition period, and carries out the smoothing by using a predetermined number of stored val- 
ues of the output data from the adaptive parameter-calculating means. 

[001 3] More preferably, the smoothing comprises calculating moving average values of the output data from the adap- 
tive parameter-calculating means. 

[001 4] Further preferably, the adaptive parameter-calculating means includes storage means for sequentially storing 
values of each of the adaptive parameters with the first predetermined repetition period, and updating means for 
sequentially updating each of the values of the each of the adaptive parameters stored in the storage means to an iden- 
tical value with the second predetermined repetition period, the smoothing means calculating moving average values 
of the values of the each of the adaptive parameters sequentially updated to the identical value. 
[0015] Also preferably, the smoothing means includes a ring buffer for storing the output data from the adaptive 
parameter-calculating means. 

[0016] Preferably, the second predetermined repetition period is set to q times the first predetermined repetition 
period (q is an integer equal to or larger than 2), the predetermined number being equal to or larger than the q. 
[001 7] Further preferably, the q depends upon a particular frequency. 

[0018] More preferably, assuming that a number of the cylinders is j. the second predetermined repetition period is 
set to j times the first predetermined repetition period, the predetermined number being equal to t times the number j of 
the cylinders (t is an integer equal to or larger than 1). 

[0019] The above and other objects, features, and advantages of the invention will become more apparent from the 
following detailed description taken in conjunction with the accompanying drawings. 

BRIEF DESCRIPTION OF THE DRAWINGS 

[0020] 

Fig. 1 is a block diagram showing the arrangement of a multi-cylinder internal combustion engine and an air-fuel 
ratio control system therefor, according to an embodiment of the invention; 

Fig. 2 is a block diagram useful in explaining a manner of controlling the air-fuel ratb of a mixture supplied to the 
engine appearing in Fig. 1 ; 

Fig. 3 is a flowchart showing a main routine for calculating an adaptive control correction coefficient KSTR in 
response to an output from a LAF sensor appearing in Fig. 1 ; 

Fig. 4 is a flowchart showing a subroutine for carrying out a LAF feedback control region-determining process, 
which is executed at a step S6 in Fig. 3; 
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Fig. 5 is a block diagram uSHPin explaining a manner of calculating the adaptive <SBc! correction coefficient 
KSTR; 

Figs. 6A and 6B are diagrams useful in explaining improved control performance, exhibited by a control method 
employed in the embodiment; 

5 Fig. 7A is a block diagram showing an adaptive parameter vector-calculating section for calculating an adaptive 
parameter vector; 

Fig. 7B is a block diagram showing a variation of the adaptive parameter vector-calculating section; 
Fig. 8A and 8B are diagrams useful for explaining characteristics of a filter appearing in Fig. 7A; 
Fig. 8C is a diagram useful for explaining characteristics of a filter appearing in Fig. 7B; 
10 Fig. 9 is a flowchart showing a subroutine for calculating the adaptive control correction coefficient KSTR; 
Fig. 10 is a flowchart showing a subroutine for calculating adaptive parameters; 

Fig. 1 1 is a diagram which is useful in explaining a method of calculating moving average values of the adaptive 
parameters; and 

Fig. 12 is a diagram showing changes in adaptive parameter values and their moving average values. 



P ENTAILED PESCRIPTION 



[0021] The invention will now be described in detail with reference to the drawings showing an embodiment thereof. 
[0022] Fig. 1 schematically shows the whole arrangement of a multi-cylinder internal combustion engine (hereinafter 

20 simply referred to as the engine) and a control system therefor, accorcfing to an embodiment of the invention. 

[0023] In Fig. 1 , reference numeral 1. .designates a four-cylinder type internal combustion engine. The engine 1 has 
an intake pipe 2 having a manifold part (intake manifold) 1 1 directly connected to the combustion chamber of each cyl- 
inder. A throttle valve 3 is arranged in the intake pipe 2 at a location upstream of the manifold part 1 1 . A throttle valve 
opening (8TH) sensor 4 is connected to the throttle valve 3. for generating an electric signal indicative of the sensed 

25 throttle valve opening 8TH and supplying the same to an electronic control unit (hereinafter referred to as the ECU) 5. 
The intake pipe 2 is provided with an auxiliary air passage 8 bypassing the throttle valve 3, and an auxiliary air amount 
control valve (electromagnetic valve) 7 is arranged across the auxiliary air passage 6. The auxiliary air amount control 
valve 7 is electrically connected to the ECU 5 to have an amount of opening thereof controlled by a signal therefrom. 
An intake air temperature (TA) sensor 8 is inserted into the intake pipe 2 at a location upstream of the throttle valve 3, 

30 for supplying an electric signal indicative of the sensed intake air temperature TA to the ECU 5. The intake pipe 2 has 
a swelled portion 9 as a chamber interposed between the throttle valve 3 and the intake manifold 1 1 . An intake pipe 
absolute pressure (PBA) sensor 10 is arranged in the chamber 9, for supplying a signal indicative of the sensed intake 
pipe absolute pressure PBA to the ECU 5. An engine coolant temperature (TW) sensor 13. which may be formed of a 
thermistor or the Ifte, is mounted in the cylinder block of the engine 1 filled with an engine coolant, for supplying an elec- 

35 trie signal indicative of the sensed engine coolant temperature TW to the ECU 5. A crank angle position sensor 14 for 
detecting the rotational angle of a crankshaft, not shown, of the engine 1 is electrically connected to the ECU 5 for sup- 
plying an electric signal indicative of the sensed rotational angle of the crankshaft to the ECU 5. The crank angle posi- 
tion sensor 14 is comprised of a cylinder-discriminating sensor, a TDC sensor, and a CRK sensor. The cylinder- 
discriminating sensor generates a signal pulse (hereinafter referred to as a CYL signal pulse) at a predetermined crank 

40 angle of a particular cylinder of the engine 1 , the TDC sensor generates a signal pulse at each of predetermined crank 
angles (e.g. whenever the crankshaft rotates through 180 degrees when the engine is of the 4<ylinder type) which each 
correspond to a predetermined crank angle before a top dead point (TDC) of each cylinder corresponding to the start 
of the suction stroke of the cylinder, and the CRK sensor generates a signal pulse at each of predetermined crank 
angles (e.g. whenever the crankshaft rotates through 30 degrees) with a predetermined repetition period shorter than 

45 the repetition period of TDC signal pulses. The CYL signal pulse. TDC signal pulse, and CRK signal pulse are supplied 
to the ECU 5, which are used for controlling various kinds of timing, such as fuel injection timing and ignition timing, and 
for detecting the engine rotational speed NE. 

[0024] Fuel injection valves 12 are inserted into the intake manifold 1 1 for respective cylinders at locations slightly 
upstream of intake valves, not shown. The fuel injection valves 12 are connected to a fuel pump, not shown, and elec- 
50 trically connected to the ECU JS to have the fuel injection timing and fuel injection periods (valve opening periods) 
th*" eof controlled by signals therefrom. Spark plugs, not shown, of the engine 1 are also electrically connected to the 
Ev. j 5 to have the ignition timing 6IG thereof controlled by signals therefrom. 

[0025] An exhaust pipe 1 6 of the engine has a manifold part (exhaust manifold) 1 5 directly connected to the combus- 
tion chambers of the cylinders of the engine 1 . A linear output air-fuel ratio sensor (hereinafter referred to as the LAF 
s$ sensor) 17 is arranged in a confluent portion of the exhaust pipe 16 at a location immediately downstream of the 
exhaust manifold 15. Further, a first three-way catalyst (immediate downstream three-way catalyst) 19 and a second 
three-way catalyst (bed-downstream three-way catalyst) 20 are arranged in the confluent portion of the exhaust pipe 1 6 
at locations downstream of the LAF sensor 1 7, for purifying noxious components present in exhaust gases, such as HC, 
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CO, and NOx. An oxygen co^Jpition sensor (hereinafter referred to as the 02 sensor^PI inserted into the exhaust 
pipe 16 at a location intermediate between the three-way catalysts 1 9 and 20. 

[0026] The LAF sensor 1 7 is electrically connected via a low-pass filter 22 to the ECU 5, fa supplying the ECU 5 with 
an electric signal substantially proportional in value to the concentration of oxygen present in exhaust gases from the 
engine (i.e. the air-fuel ratio). The 02 sensor 18 has an output characteristic that output voltage thereof drastically 
changes when the air-fuel ratio of exhaust gases from the engine changes across a stoichiometric air-fuel ratio to 
deliver a high level signal when the mixture is richer than the stoichiometric air-fuel ratio, and a low level signal when 
the mixture is leaner than the same. The 02 sensor 18 is electrically connected via a low-pass filter 23 to the ECU 5 for 
supplying the ECU 5 with the high or low level signal. The low-pass filters 22 and 23 are provided for eliminating high 
frequency noise components, and influence thereol on the responsiveness of the air-fuel ratio control system is negli- 
gible. 

[0027] The engine 1 includes a valve timing changeover mechanism 60 which changes valve timing of at least the 
intake valves out of the intake valves and exhaust valves, not shown, between a high speed valve timing suitable for 
operation of the engine in a high speed operating region thereof and a low speed valve timing suitable for operation of 
is the engine in a low speed operating region thereof. The changeover of the valve timing includes not only timing of open- 
ing and closing of the valve but also changeover of the valve lift amount, and further, when the low speed valve timing 
is selected, one of the two intake valves is disabled, thereby ensuring stable combustion even when the air-fuel ratio of 
the mixture is controlled to a leaner value than the stoichiometric air-fuel ratio. 

[0028] The valve timing changeover mechanism 60 changes the valve timing by means of hydraulic pressure, and an 
20 electromagnetic valve for changing the hydraulic pressure and a hydraulic pressure sensor, neither of which is shown, 
are electrically connected to the ECU 5. A signal indicative of the sensed hydraulic pressure is supplied to the ECU 5 
which in turn controls the electromagnetic valve for changing the valve timing. 

[0029] Further electrically connected to the ECU 5 is an atmospheric pressure (PA) sensor 21 , for detecting atmos- 
pheric pressure PA. and supplying a signal indicative of the sensed atmospheric pressure PA to the ECU 5. 

25 [0030] The ECU 5 is comprised of an input circuit having the functions of shaping the waveforms of input signals from 
various sensors including ones mentioned above, shifting the voltage levels of sensor output signals to a predetermined 
level, converting analog signals from analog-output sensors to digital signals, and so forth, a central processing unit 
(hereinafter referred to as the CPU), a memory circuit comprised of a ROM storing various operational programs which 
are executed by the CPU and various maps and tables, referred to hereinafter, and a RAM tor storing results of calcu- 

30 lations from the CPU, etc., and an output circuit which outputs driving signals to the fuel injection valves 12 and other 
electromagnetic valves, spark plugs, etc. 

[0031 ] The ECU 5 operates in response to the above-mentioned signals from the sensors to determine operating con- 
ditions in which the engine 1 is operating, such as an air-fuel ratio feedback control region in which air-fuel ratio feed- 
back control is carried out in response to outputs from the LAF sensor 1 7 and the 02 sensor 18. and air-fuel ratio open- 
35 loop control regions, and calculates, based upon the determined engine operating conditions, the fuel injection period 
TOUT over which the fuel injection valves 12 are to be opened, by using the following equation (1 ), to output signals for 
driving the fuel injection valves 12, based on results of the calculation: 



40 



TOUT = TIMF x KTOTAL x KCMDM x KFB (1) 



where TIMF represents a basic value of the fuel injection amount TOUT, KTOTAL a correction coefficient, KCMDM a 
final desired air-fuel ratio coeff icient, and KFB a feedback correction coefficient, respectively. 
[0032] Fig. 2 is a block diagram useful in explaining the manner of calculating the fuel injection period TOUT by using 
the equation (1). With reference to the figure, an outline of the manner of calculating the fuel injection period TOUT 
45 according to the present embodiment will be described. It should be noted that in the present embodiment, the amount 
of fuel to be supplied to the engine is calculated, actually, in terms of a time period over which the fuel injection valve 6 
is opened (fuel injection period), but in the present specification, the fuel injection period TOUT is referred to as the fuel 
irqection amount or the fuel amount since the fuel injection period is equivalent to the amount of fuel injected or to be 
injected. 

so [0033] In Fig. 2, a block 81 calculates the basic fuel amount (basic value of the fuel injection amount) TIMF corre- 
sponding to an amount of intake air supplied to the engine 1 . The basic fuel amount TIMF is basically set according to 
the engine rotational speed NE and the intake pipe absolute pressure PBA. However, it is preferred that a model repre- 
sentative of a part of the intake system extending from the throttle valve 3 to the combustion chambers of the engine 1 
is prepared in advance, and a correction is made to the basic fuel amount TIMF in dependence on a delay of the flow 

55 of intake air obtained based on the model. In this preferred method, the throttle valve opening GTH and the atmospheric - 
pressure PA are also used as additional parameters indicative of operating conditions of the engine. 
[0034] Reference numerals B2 to B4 designate multiplying blocks, which multiply the basic fuel amount TIMF by 
respective parameter values input thereto, and deliver the product values. These blocks carry out the arithmetic opera- 
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8'Wtb 0 (k).r 1 (k)....r nH<M (k).s 0 (k) Sn .,(k)] ^ (4) 

=lb 0 (k).r 1 (k).i 2 (k).f 3 (k),8 0 (k)] 

C T (k)=[u(k) u(k-n*d+1). y(k) y(k-n+1)] (5) 

=[u(k). u(k-1), u(k-2). u(k-3). y(k)] 

[0052] The adaptive parameter vector e(k) is expressed by the following equation (6): 

e(k)=6(k-l)+r(k-l)C(k-d)e'(k) (6) 

where the symbols r(k) and e*(k) represent a gain matrix and an identification error signal, respectively, and can be 
expressed by the following recurrence formulas (7) and (8): 

m . . r M . , W WK^r(K.i) 1 (?) 



- -e-(k)= (D ( z1) yW' 6 V 1 )^M) (8) 
1+C T (k-d)r(k-1)C(k-d) 

[0053] Further, rt is possible to provide various specific algorithms depending upon set values of X t (k) and in the 
equation (7). For example, if Xj (k) = 1 and X 2 (k) «= X (0 < X < 2) hold, a progressively decreasing gain algorithm is pro- 
vided (if X = 1 , the least square method), if X , (k) = X , (0 < X 1 < 1 ) and X 2(k) = X 2 (0 < X^ < 2) hold, a variable gain 
algorithm Of = 1 t the method of weighted least squares), and if x 1 (k)/X 2 (k) = a and if X3 is expressed by the follow- 
ing equation (9), X ^k) =: X 3 provides a fixed trace algorithm Further, if X,(k) = 1 and Xg{k) = 0 hold, a fixed gain algo- 
rithm is obtained. In this case, as is dear from the equation (6). r(k) = r(k-1) holds, and hence r(k) = r (fixed value) 
is obtained. 

[0054] Further. D(Z" 1 ) in the equation (8) is an asymptotically stable polynomial which can be defined by a system 
designer as desired to determine the convergence of the system. In the present embodiment, it is set to a value 1 .0. 

a+C T (kKi)r(k-1)C(k-d) trr <°> 

[0055] In the equation (9), trF(0) is a trace function of the matrix r(0). and specifically, it is a sum (scalar) of diagonal 
components of the matrix r(0). 

[0056] The STR controller and the adaptive parameter-adjusting mechanism are arranged outside the fuel injection 
amount-calculating system, and operate to calculate the adaptive control correction coefficient KSTR(k) such that the 
actual equivalent ratio KACT(k+d) becomes equal to the desired equivalent ratio KCMD(k) in an adaptive manner. 
[0057] As shown in Fig. 5, the adaptive control correction coefficient KSTR(k) and the actual equivalent ratio KACT(k) 
are determined, which are input to the adaptive parameter-adjusting mechanism, where the adaptive parameter vector 
6(k) is calculated to be input to the STR controller. The STR controller is also supplied with the desired equivalent ratio 
KCMD(k) and calculates the adaptive control correction coefficient KSTR(k) such that the actual equivalent ratio 
KACT(k+d) becomes equal to the desired equivalent ratio KCMD(k). by using the following recurrence formula (10): 

ueTD . w KCMD(k)-s 0 xKACT(k)-r ^KSTRfk-IJ-r 2 xKSTR(k-2)-r 3 xKSTR(k-3) 

K5TR(k)= - (10) 

D 0 



[0058] The equation (10) (and the equations (4) and (5)) is obtained when the STR controller is designed on assump- 
tion that the dead time d of the engine 1 as a plant and the LAF sensor 17 as means for detecting the output from the 
plan is as long as three control cycles. However, due to a modification of the specifications of the engine 1 or the LAF 
sensor 17. in some cases, the dead time d can be longer than three control cycles. Assuming that the dead time d = 5 
holds, for instance, the adaptive parameter vector 8(k) and the input £T(k) to the adaptive parameter-adjusting mecha- 
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nism are expressed by the 



equations (11) and (12): 



5 



e T (k)=[b 0 (k), r t (k), r 2 (k), r 3 (k). r 4 (k), r 5 (k). s 0 (k)] 
C T (k)=[u(k), u(k-1), u(k-2), u(k-3), u(k-4), u(k-5). y(k)] 



(11) 



(12) 



[0059] Further, the equation for calculating the adaptive control correction coefficient KSTR(k) is expressed by the 
following equation (13): 



is [0060] Therefore, the amount of arithmetic operations required to calculate the adaptive control correction coefficient 
KSTR largely increases to such a level as will make it impractical to calculate the coefficient KSTR by using the CPU 
installed on the automotive vehicle for control of the engine. To eliminate this inconvenience, in the present embodi- 
ment, the configuration of the STR controller is adapted to the dead time d = 3 as shown in Fig. 5, and sampling timing 
(in the present embodiment, the operation cycle in which is calculated the input vector < used in the arithmetic opera- 

20 tions is referred to as sampling timing) of the input vector X, to the adaptive parameter-adjusting mechanism is set to a 
period corresponding to an actual dead time DACT (DACT > d, e.g. 5) (hereinafter referred to as "lowered-order STR"). 
More specifically, in the above equations (6) to (8). if d = DACT = 3 holds . £(k-d) is made equal to C(k-3). which is also 
the case with t T (k-d). i.e. the transposed determinant of t(k-d), and the same applies to the following description. How- 
ever, if the actual dead time DACT is larger than 3, e.g. equal to 5, the STR controller per se is adapted to d = 3, and 

25 the input vector <(k-d) of the equations (6) to (8) is set to £(k-5) which is adapted to the actual dead time DACT. Then, 
the adaptive parameter vector 6(k) (which is composed of elements bo. Sq. r 1§ r 2 . and r 3 exclusive of r 4 and r s ) is applied 
to the equation (10) to calculate the adaptive control correction coefficient KSTR. In this case, the vector t;(k) is 
expressed by the the right side of the equation (5). which naturally corresponds to d = 3, and the orders of terms of the 
equation applied are not increased. The employment of the lowered-order STR makes it possible to carry out the adap- 

30 tive control in a manner adapted to the actual dead time DACT without increasing the orders of terms of equations used 
in the STR controller and. the adaptive parameter-adjusting mechanism, and hence carry out high-accuracy adaptive 
control while minimizing an increase in the amount of arithmetic operations required to obtain the adaptive control cor- 
rection coefficient KSTR. 

[0061 ] Figs. 6A and 68 show examples of changes in the adaptive correction coefficient KSTR and the actual equiv- 
35 alent ration KACT which occurred when the desired equivalent ratio KCMD was changed in the case where the actual 
dead time DACT of the controlled object is longer than three control cycles. Fig. 6A shows a case in which the STR con- 
troller is adapted to the dead time d = 3, and the sampling timing of the input vector £T(k) to the adaptive parameter- 
adjusting mechanism is also adapted to the dead time d = 3 by setting the factor £(k-d) in the equations (6) to (8) to C(k- 
3), while Fig. 6B shows a case in which the STR controller is adapted to the dead time d = 3, and the sampling timing 
40 of the input vector C to the adaptive parameter-adjusting mechanism is adapted to the dead time d«4, which is sub- 
stantially equal to the actual dead time DACT by setting the factor £(k-d) in the equations (6) to (8) to £(k-4). As is clear 
from these figures, even if the dead time d employed in designing the STR controller is shorter than the actual dead 
time DACT, it is possible to largely improve the responsiveness and stabifity of the adaptive control. 
[0062] In general, the actual dead time of the controlled object is an analog value which changes continuously. How- 
45 ever, in the adaptive control, the dead time has to be converted into a discrete or digital value. Therefore, assuming that 
the actual dead time is as long as 4.5 control cycles, for instance, it is preferable to select the more suitable one of val- 
ues 4 and 5 for the dead time through comparison between control performance exhibited when the dead time DACT 
is set to 4 and that exhibited when the same is set to 5. 

[0063] As described above, the employment of the lowered-order STR can enhance the performance of the adaptive 
so control without increasing the amount of arithmetic operations. However, the model assumed in designing the control 
system is different from the actual object to be controlled, so that the adaptive parameter vector 6 can be liable to drift. 
Further, even if the lowering of orders of equations is not carried out. the equation (6) performs accumulation of slight 
identification errors caused by external disturbances in the adaptive parameter vector. Therefore, so long as the adap- 
tive parameter vector 0 is calculated by using the equation (6). it is inevitable that the adaptive parameter vector e drifts 
55 when the engine continues to be in a steady operating condition. 

[0064] To eliminate this inconvenience, in the present embodiment, means is additionally provided for preventing the 

adaptive parameter vector from drifting. This means wil be described in detail hereinafter. 

[0065] Fig. 7A shows an adaptive parameter vector-calculating section employed in the control system according to 
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al^Wg the adaptive parameter vector 9. In the figure, ei^^Hr* 



the present embodiment for cal^Pig the adaptive parameter vector 9. In the figure, eid|Pfnd Q(k) represent the 
identification error and the variable gain, respectively, which are defined by the following equations (14) and (15), 
respectively: 

eid(k)=D(Z')y(k)-e\k.m*-<t) (14) 

OH. r ( k.i)^) 
i+C T (k-d)r(k-i)C(k^j) 

[0066] The identification error eid(k) corresponds to the numerator of the identification error signal e* (equation (8)), 
while the variable gain Q(k) corresponds to a quotient obtained by dividing the second term of the equation (6) by the 
identification error eid(k). When the identification error eid and the variable gain Q(k) are used, the adaptive parameter 
vector e(k) is expressed by the following equation (16): 

6(k)=8(k-1)+Q(k)eid(k) (16) 

[0067] As shown in Fig. 7A, the identification error eid(k) is input to a nonlinear filter B51 and an output eida(k) there- 
from is input via a multiplier B52 to an integrator B53. The nonlinear filter B51 has an input/output characteristic e.g. as 
shown in Fig. 8 A. More specif icafly, when -r\ ss eid £ r\ holds (n represents a predetermined value empirically obtained), 
which means that eid falls within the deadzone defined by -tj and tj , the output eida is set equal to 0 and when the iden- 
tification error eid is outside the dead zone, the output eida is set equal to the input eida. 

[0068] The arithmetic operations carried out at the blocks B52 and B53 correspond to those carried out by using the 
equation (16) in which eid(k) is replaced by eida(k). 

[0069] In the present embodiment, the process of the nonlinear filter B51 is additionally carried out. As a result, when 
the engine enters a steady operating condition, and identification of the adaptive parameter vector 6 is substantially 
completed, the identification error eid falls within the dead zone defined between -n and +n. so that the output eida from 
the block B51 is equal to 0, which prevents slight identification errors caused by differences in characteristics between 
the model assumed in designing the control system and the actual plant and external disturbances from being accumu- 
lated in the adaptive parameter vector 9, to thereby prevent drifting of the adaptive parameter vector 6. 
[0070] The input/output characteristic of the nonlinear filter B51 shown in Fig. 8A shows a discontinuity as the input 
eid changes across -q, or tj, which causes a drastic change in the vector 0 when the identification error eid goes beyond 
the dead zone. To eliminate this inconvenience, it is preferable that the nonlinear filter B51 has an input/output charac- 
teristic free from Discontinuities as shown e.g. in Hg. 8B. That is, when -r\ * eid ^ x\ holds, the output eida is set to 0. 
while when eid > i\ holds, the same is set to eid -r\, and when eid < -rj, holds, the same is set to eid + r\. 
[0071] Further, the arrangement of Rg 7A in which the nonlinear filter B51 alone is added undergoes occurrence of 
a steady state error of the adaptive control (steady state difference between the desired equivalent ratio KCMD and the 
actual equivalent ration KACT). This is because the identification error within the dead zone does not reflect on the 
adaptive parameter vector 6. 

[0072] To overcome this inconvenience, it is desirable that blocks B61 to B66 are additionally provided, as shown in 
Fig. 7B. Multiplies B61, B64, an adder B62, and a delay B63 constitute a first-order lag filter (low-pass filter). The output 
eidf from this filter has Hs level limited by a limiting filter B65, and is added to the output eida from the nonlinear filter 
B51 by an adder B66. The output from the adder B66 is input to the multiplier B52. 
[0073] The output eidf from the first-order lag filter eidf is expressed by the following equation (17): 

eidf = CLF x eidfa(k-l) + (1 - CLF) x eid(k) (1 7) 

wherein CLF represents an averaging coefficient which is set to a value smaller than 1 but very close to 1, e.g. "0.988*. 
In other words, the coefficient CLF is provided to set the cutoff frequency of the low-pass filter to a value close to 0. 
[0074] The limiting filter 65 constitutes a limiter having a characteristic as shown in Rg. 8C. More specifically, when 
the output eidf from the first-order lag filter, which is input to this limiter, is smaller than -r\, an output eidfa from the limiter 
is set to -tj. whereas when the output eidf is larger than r\. the output eidfa is set to r\. When -tj * eid * r\ holds, eidfa is 
set to eidf. 

[0075] In a variation shown in FIG. 7B, a steady state error component of the identification error eid is extracted by 
the low-pass fitter having a cut-off frequency dose to 0 and added to the output eida from the nonlinear filter B51 . As a 
result, the steady state error component in the identification error eid is reflected on the adaptive parameter vector 6 to 
thereby reduce the steady state error of the adaptive control (steady state error between the desired equivalent ratio 
KCMD and the actual equivalent ratio KACT). 
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[0076] Further, the provisicjjM^e limiting filter B65 prevents the steady state error ^Jbnent of the identification 
error eid used for calculation of the adaptive parameter vector 0 from becoming larger than the original value due to the 
addition by the adder B66 when the identification error falls outside the dead zone (i.e. eid > 11 or eid < -x\). 
[0077] Next the equation for calculating the adaptive control correction coefficient KSTR actually employed in the 
present embodiment will be descrtoed. The above equations (5) to (10) are applied to a case where the control cycle 
and the repetition period of calculation of the KSTR value (repetition period of generation of TDC signal pulses) coincide 
with each other and the adaptive control correction coefficient KSTR thus calculated is commonly used for all the cylin- 
ders. In the present embodiment, however, the control cycle is as long as four TDC signal pulses corresponding to the 
number of cylinders, whereby the adaptive control correction coefficient KSTR is determined cylinder by cylinder. More 
specifically, the above-mentioned equations (5) to (10) are replaced by the following equations (18) to (23), respectively, 
to calculate the adaptive control correction coefficient KSTR cylinder by cylinder for use in the adaptive control: 

C T (k)=[u(k), u(k-4), u(k^8), u(k-12), y(k)] (18) 

6(k)=e(k-4)+r(k-4)<(k-4xd)e*M ( 19 > 

1 T X 2 (k)r(k-4)C(k-4xd)C T (k-4xd)r(k-4)1 

Hk)=r7d r(k-4)- 21 ' " V T — — (20) 

M k) L X 1 (k)+X 2 (k)C T (k-4xd)r(k-4)C(k-4xd)J 

^ Pb'MyW-e^tMxd) (21) 
1+t T (k-4xd) r(k-4)C(Mxd) 



2 



Mk)=1 ir(k-4)^k-4xd)r . i m 

a+C T (k-4xd)r(k-4)ftk-4xd) trr (0) 

KSTR(k)= {KCMD (k)-s 0 x KACT (k)-r \ x KSTR (k-4) -r 2 x KSTR (k-8)-r 3 x KSTR (k-12)}/b 0 (23) 

[0078] It should be noted that when the actual dead time DACT is e.g. 4, the adaptive parameters bo. Sq, to r 3 are 
calculated by using the equations (19) to (22) with d set to 4. 

[0079] Fig. 9 shows a subroutine for calculating the adaptive control correction coefficient KSTR which is executed at 
the step S9 in Fig. 3. 

[0080] First, at a step S401 , it is determined whether or not the reset flag FKLAFRESET assumed "1" in the last loop. 
If FKLAFRESET = T held in the last loop, which means that the adaptive control was not carried out in the last loop, 
the initialization of the adaptive parameters bo, Sq and r t to r 3 is carried out by setting them to respective initial values, 
and then the program proceeds to a step S404. On the other hand, if it is determined at the step S401 that FKLAFRE- 
SET = "1 " held, which means that the adaptive control was carried out in the last loop as well, the process for calculating 
the adaptive parameters b 0 . Sq and r, to r 3 is carried out by executing a routine shown in Fig. 10. 
[0081] In the present embodiment, the calculation of 6 (k), i.e. adaptive parameters bo. Sq and r n to r 3 by using the 
equation (19) is carried out once per four TDC periods (time period over which four TDC signal pulses are generated, 
i.e. one combustion cyde).Therefore. at a step S431 in Fig. 10. it is determined whether or not four TDC periods have 
elapsed from the last calculation of the adaptive parameters using the equation (19). If it is determined that four TDC 
periods have elapsed, the adaptive parameters t>o(k) t so(k). and r t (k) to r 3 (k) are calculated by using the equation (19) 
at a step S432. If four TDC periods have not elapsed, the adaptive parameters bo. Sq and r\ to r 3 are set to respective 
immediately preceding values bo(k-1), So(k-1). and r^k-1) to r 3 (k-1). 

[0082] After execution of the step S432 or S433, moving average values boAV, SqAV, r t AV, r 2 AV, and r 3 AV over p TDC 
periods (e.g. p = 8, i.e. 8 TDC periods) are calculated by using the following equations (24) to (28) at a step S434, fol- 
lowed by terminating the program. The control system includes a ring buffer (memory means) for storing values of the 
adaptive parameters b 0 . Sq, to r 3 obtained over the p TDC periods (hereinafter p will be referred to as "the averaging 
perkxT) to calculate the moving average values. The contents of the ring buffer are updated to newly-calculated adap- 
tive parameter values or the immediately preceding values whenever one TDC period elapses, whereby the oldest data 
stored therein are erased. The updating of the contents of the ring buffer may be carried out by storing a newly calcu- 
lated value of the adaptive parameter vector 6 over four TDC periods (i.e. by storing the newly calculated value e(k) as 
the values 6(k) to 6(k+3): 
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connected to said cylinders^BRjing air-fuel ratio-detecting means arranged in said SBBst system for detecting 
an air-fuel ratio of a mixture supplied to said engine and for generating an output indicative of said air-fuel ratio of 
said mixture, an adaptive controller for determining an amount of fuel to be supplied to said engine with a first pre- 
determined repetition period in a manner such that said output from said air-fuel ratio-detecting means becomes 
equal to a desired value, and adaptive parameter-adjusting means for adjusting adaptive parameters used by said 
adaptive controller, 

the improvement wherein: 



said adaptive parameter-adjusting means comprises: 

adaptive parameter-calculating means for calculating said adaptive parameters with a second predetermined 
repetition period longer than said first predetermined repetition period and tor generating output data indicative 
of results of said calculation; and 

smoothing means for smoothing said output data from said adaptive parameter-calculating means and for gen- 
erating output data indicative of the smoothed data with a repetition period at least equal to said first predeter- 
mined repetition period; 

wherein said adaptive controller uses the smoothed data generated from said smoothing means as values of 
said adaptive parameters. 

2. An air-fuel ratio control system according to claim 1 , wherein said smoothing means stores said output data from 
said adaptive parameter-calculating means with said first predetermined repetition period, and carries out said 
smoothing by using a predetermined number of stored values of said output data from said adaptive parameter- 
calculating means. 

3. An air-fuel ratio control system according to claim 2. wherein said smoothing comprises calculating moving average 
values of said output data from said adaptive parameter-calculating means. 

4. An air-fuel ratio control system according to claim 3. wherein said adaptive parameter-calculating means includes 
storage means for sequentially storing values of each of said adaptive parameters with said first predetermined 
repetition period, and updating means for sequentially updating each of said values of said each of said adaptive 
parameters stored in said storage means to an identical value with said second predetermined repetition period, 
said smoothing means .calculating moving average values of said values of said each of said adaptive parameters 
sequentially updated to said identical value. 

5. An air-fuel ratio control system according to claim 2, wherein said smoothing means includes a ring buffer for stor- 
ing said output data from said adaptive parameter-calculating means. 

6. An air-fuel ratio control system according to claim 2, wherein said second predetermined repetition period is set to 
q times said first predetermined repetition period (q is an integer equal to or larger than 2), said predetermined 
number being equal to or larger than said q. 

7. An air-fuel ratio control system according to daim 6, wherein said q depends upon a particular frequency. 

8. An air-fuel ratio control system according to claim 3, wherein said second predetermined repetition period is set to 
q times said first predetermined repetition period (q is an integer equal to or larger than 2), said predetermined 
number being equal to or larger than said q. 



9. An air-fuel ratio control system according to daim 2, wherein assuming that a number of said cylinders is j. said 
second predetermined repetition period is set to j times said first predetermined repetition period, said predeter- 
mined number being equal to t times said number j of said cylinders (t is an integer equal to or larger than 1 ) . 

10. An air-fuel ratio control system according to daim 3, wherein assuming that a number of said cylinders is j. said 
second predetermined repetition period is set to j times said first predetermined repetition period, said predeter- 
mined number being equal to t times said number j of sad cylinders (t is an integer equal to or larger than 1). 
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